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Abstract 

The estimation of rare K decay matrix-elements from Ki^ experimental data is extended beyond LO in 
Chiral Perturbation Theory. Isospin-breaking effects at NLO (and partially NNLO) in the ChPT expansion, 
as well as QED radiative corrections, are now accounted for. The analysis relies mainly on the cleanness of 
two specific ratios of form-factors, for which the theoretical control is excellent. As a result, the uncertainties 
on the — > Tr^iyp and Kl — > tt'^i^P matrix-elements are reduced by a factor of about 7 and 4, respectively, 
and similarly for the direct CP-violating contributions to — > Tr'^e^e" and Kl —> Tr^/i^/^". They could 
be reduced even further with better experimental data for the Ki^ slopes and the K^^ branching ratios. 
As a result, the non-parametric errors for B (K — s- tti^P) and for the direct CP-violating contributions to 
B [Kl — > 7r"^+£^) are now completely dominated by those on the short-distance physics. 



1 Introduction 



The rare decays K nvP and Kl 'K^t^t' are important tools to test the Standard Model (SM), and 
to search for possible new physics. As they proceed through flavor changing neutral currents (FCNC), 
though they are very suppressed in the SM, they show an exceptional sensitivity to short-distance physics. 
Recent theoretical progress has greatly improve the theoretical control over the SM predictions [1-3]. On 
the experimental side, rare K programs at CERN and J-PARC are currently under study, aiming at a 
measurement by the beginning of the next decade. 

An essential aspect of these decay modes is of being semi-leptonic. Contrary to e'/e for example, the dom- 
inant contribution does not come from four-quark operators, but rather from quark vector-current FCNC 
operators like {s^^l,>}){y^^ {\ — ^-^v), on which we have an excellent control. This stems from their rela- 
tionship under the isospin symmetry with the charge-current (CC) Fermi operators {s'y^u){i'y^ (1 — 75) i/e). 
More specifically, the FCNC hadronic matrix-elements required for rare K decays as well as the one for the 
CC transition in K ^ wlu^ decays {K^^) are parametrized as {T = P — K) 

{tt^ [K] \qr\aq\K\P)) = -^ (T^) [P + Kf + JK^-^T^) {P ^ Kf ) , (1) 

with dj some Clebsch-Gordan coefficients and the Gell-Mann matrices Xfcnc = Ae ± iXr, Xcc = A4 ± iAs 
projecting out the desired quark-flavor structures. In the isospin limit, all these form-factors are equal. 

At lowest order in Chiral Perturbation Theory (ChPT), these form- factors arc all related to the conserved 
current form-factors f^''~ {T"^), and thus ff''\T^) = 1, f^'^\T^) = for all i,j = -h,0. Further, the 
Ademollo-Gatto theorem protects against large SU (3) corrections, which can arise only at the second order 
in nis — niu^d- Of course, in practice, it is not useful for us to compute these SU (3) corrections since they are 
the same for all K ^ n transitions. Indeed, our goal is to use instead the precise experimental information 
on the form- factors obtained from K^^ decays. The problem then reduces to the study of isospin-breaking 
effects, proportional to e^^^ ~ (m„ — md)/ms, to relate FCNC and CC form-factors as precisely as possible. 

This strategy of using K^s data is common practice, but currently relies on the 0(p^£*^^)), LO analysis 
of Ref. [4]. Given the recent theoretical progress in the computation of short-distance QCD corrections, it 
is now time to improve and go beyond LO. More precisely, with respect to Ref. [4], our goals are: 

1 - To include isospin-breaking effects at NLO (and partially NNLO) in the ChPT expansion. 

2 - To account for QED radiative corrections, at leading order in the ChPT expansion. 

3 - To update matrix-elements using the latest K^^ experimental data. 

4 - To perform a detailed error study, including both theoretical and experimental uncertainties. 

The outline of the paper is as follows. In the next Section, the master formulas arc given. In Section 
3, the ChPT results for the form-factors at 0{p^e'^^^), and partially 0{p^s'^^^), are given and discussed. 
Special emphasis will be set on two ratios of form-factors on which an exceptional theoretical control can 
be reached. Then, the numerical analysis is performed in Section 4 and our results are summarized in the 
Conclusion. Finally, loop functions as well as the details of the computation of the QED radiative corrections 
are presented in the Appendix. 

2 Generalities 

The FCNC weak Hamiltonian relevant for the K ttvv decays and for the direct CP-violating contribution 
to the Kl — > 7r°^+^~ decays is 

Heff = )= — - > ^—r^T-Qt.u 1= > {yrvQrv + ViaQia) + h-c. , (2) 

Q.U = {sj^^d) X {va" (1 - 75) yt) , Qiv = {sj^^d) x (i^'^e) , QrA = (s7^rf) x [h'^j^i) ■ (3) 
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In the SM, the Wilson coefScient y^, is given by 

= (ReAt + ilmXt) Xt + \Vusf ReAcP„,c , (4) 

with Xg = Xt = 1.464 ± 0.041 [1], P„,c = 0.41 ± 0.04 for rric (rric) = 1.30 ± 0.05 [1, 2], while y7v and 

y7A are given by [5] 

yrv (m « 1 GeV) = (0.73 ± 0.04) ImAt, yrA {Mw) = (-0.68 ± 0.03) ImXt . (5) 

In the final numerical applications, we will also use the CKM parameters as obtained in Ref. [6] (compatible 
with Ref. [7]), and from Ref. [8]: 

Im At = IMSto.lll, R-e At = -Smto.H, Re Ac = -0.2209ltg;ggg^3, = 0.2245 ± 0.0016 (6) 

All the short-distance physics is encoded in the Wilson coefficients. The remaining task for computing 

the branching ratios is thus to get the matrix-elements of these operators, and to carry out the phase-space 
integration. For later use, we collect below the parametrizations of the branching ratios in terms of the 
form-factors. For that, we adopt the standard parametrization in terms of the form-factor values at the 
origin (T^ = 0) and of their derivatives there, i.e. their slopes. In this respect, remember also that f^''^' 
can be eliminated in favor of the scalar form- factors (matrix-elements of the scalar current qXaq)'- 

{T^) = ff^' {T') + jp^^f^'""' (T') , (7) 
such that only the vector form-factors at the origin are needed, fo^ '^' (0) = fl^'^^' (0). 



Decay rates 

The Kes decay rates are 



pOEW^Ki 



K«x/f-^ (0)|'jf (l + 2A^^^^) , 



(8) 



with the Clebsch-Gordan coefficients Co+ = 1 and C+o = l/\/2. The Fermi constant Gp is fixed from 

/i decay as Gp = 1.166371 (6) • 10"" McV"^ [gj Correspondingly, the factor Spw = 1.0232 stands for 
the short distance part of the corrections to the semi-leptonic weak charged current, to leading order in 
'mz,w 00 [10]. The phase-space integrals I^-' , which are functions of the form-factor slopes, are given 
below. 

Adopting a similar parametrization, the —* ir^vv and Kl ir^yy branching ratios are written as 



B {K+ n+vv (7)) = K+ (1 + Abm) 



W,, 



B Kr 



Im yi, 



GlMl,^,a{Mzf 



2567r5 sin"^ Ow 



X /r (0) 



We will further use [1] 



^Ms(^^)"' = 127.9, (sin^ ^^^)_ = 0.231 



(9) 
(10) 

(11) 



when computing The experimental errors on a-jjg (Mz) and (sin'^ ^^)~ms "^^^^ uncertainties 

related to (unknown) higher-order, purely electroweak short-distance corrections are understood to be ac- 
counted for in the Wilson coefficients [11]. Finally, the direct CP- violating contribution to the Kl — > ■jt^£~^£~ 
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decays, which arises from the Qta and Qrv operators (for the general structure and a description of the 
indirect CP- violating and CP-conserving contributions, see for example [3]), takes the form 



DCPV 



Im yyy 



IK, 



\Vus 



3847r3 



IK, 



Vus X ff-" (0) 



■ji 
±11 . 



(12) 



(13) 



In all cases, the long-distance electromagnetic corrections are moved into the A^-'^^^'s, which thus include 
both virtual photon exchanges and real photon emissions. Because of the latter, these corrections depend 
on the experimental requirements enforced on the kinematics of the corresponding radiative decay. 

For i^^s, also the local QED corrections (including the electromagnetic and semi-leptonic counterterms) 
and the QED corrections to the phase-space arc understood in A^"*^^^. These have been computed both 
restricting the ii'^s^ phase-space to the three-body kinematics [12, 13], and for the fully inclusive case [14, 15]. 

For TT+i/P, the correction factor Aem obviously concerns the hadronic part only, and will be 

computed later on. Finally, for 'iT^i~^£~, long-distance QED radiative corrections concern only the 

lepton pair and are neglected. 

Phase-space integrals 

We use the quadratic and linear parametrizations for the vector and scalar form- factors, respectively: 

f^' = f^' (0) (l + Af + Af , /o"'-^ id') = /f (0) (l + A--^) . (14) 

Then, the phase-space integrals I*^ depends only on the slopes A'^ . Explicitly, the K^s phase-space integrals 
are {i,j = +,0 or 0, +) 



= / ..(l-i) A. 



A^ 1 + 



2z 



+ 



fr^' (0) 



n 

z 



(15) 



with z = /M'^i, ri^.^ = m^^^j/M^i, A,r = A^/^ (l,^;,r^), while those for the rare K decays are 



dzXt 



fK^.' (0) 



4r? V Z 



{i = +,L), 



(Z) 



n^" (0) 



AMl-M 



/r (0) 



ff''" (z) 



(16a) 
(16b) 
(16c) 



with = y/T—4rf/z (see Appendix A for numerical expressions). 



3 Vector form-factors for K ^ n transitions 

As we will see in the next section, devoted to the numerical analysis, the accuracy of Ki^ data are now 
below the percent level. To make full use of these impressive results, the theoretical analysis has to reach 
a corresponding level of precision. This is not an easy task as it requires some control over the 0{p^) 
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corrections. Fortunately, for the purpose of estimating FCNC form-factors from those extracted from K^s 
data, two well-chosen ratios are sufficient [4]. Considering ratios has the immediate advantage that only 
isospin-brcaking corrections, proportional to e*^^^ ^ (m„ — md)/ms, have to be brought under control. 

Therefore, in a first stage, we have extended the O [p^e'-'^^^p^a) analysis of Ref. [12] to the FCNC form- 
factors. In a second stage, the local corrections of Oip^e^'^^) to the ratios will be discussed. The good surprise 
is that the structure of the K ^ w vector transitions somehow protects our two ratios from large corrections. 

3.1 The vector form-factors at O [p'^e'^'^^-.p^a) 

For the {Tr''\V^\K'') vector form-factors at O (p'^ e^^^ : aj , one has to compute meson and photon one-loop 
diagrams, with vertices coming from the leading O (p^c^^^) cliiral Lagrangian and from the 0{p'^a) Dashen 
term Ze^e^i^^ (U^QUQ') [16]. Renormalization is carried out by adding the contributions from the 0{p'^£^'^^) 
strong counterterms Lj [17] and from the 0{p'^a) QED counterterms Ki [18]. In all these Lagrangians, the 
CC and FCNC vector currents are introduced through the covariant derivative with a A4 ± lAs and Ae ± 1X7 
flavor structure, respectively, or directly through their contributions to the field-strengths. 

In our computation, we keep the vector currents as external sources, to be coupled later to the lepton 
pairs, so that the semi-leptonic counterterms [19] are not needed to got a UV-finite result. In that respect, 
the situation is different for Ki-^ and rare K decays. For the former, these counterterms are ultimately 
needed to renormalize the QED corrections where a photon is exchanged between the lepton and the meson. 
Altogether, these form the long-distance QED corrections and are included in Aem- On the other hand, for 
rare K decays, there is no such photon exchange, and no semi-leptonic counterterm is needed. Further, it 
should be clear that short-distance QED corrections are to be accounted for in the Wilson coefficients of the 
effective operators. In particular, the Sew factor does not occur for rare K decays (see also the discussion 
in Refs. [4,12]). 

A straightforward computation gives the following {Tr^\V^\K'') matrix-elements, including terms up to 

ff-° (g2) = 1 + 3ij^„ + ff-" + 2Hl\ - V3£(2) (1 + + Sff^" + 2HJ,l) - ^/3e(4) , (I7a) 

/^+-° (q^) = 1 + 3Hl^ + iJ-'; + 2iJ]^; + V3e(2)(l + + 3i7]^+ + 2i/]^:) + VSe'^) - ^^Im , (17b) 

ff-^ = 1 + 3Hl^ + H],\ + 2HlX + 2v/3£(2)(<+ - Hl^) - ^S^m , (17c) 

ff-^ = 1 + SHI, + H^l + 2Hl\ - 2x/3e(2)(if^o - Hi,) - ^Jem {q\ M%^,M^^) , (17d) 
with 

ShM = l- + fog ^ - I log ^ . (18) 

The loop functions are given in the appendix. The two-point functions H^ involve Lg, while the three-point 
function Jem is free of counterterms but IR-divergent. The isospin-breaking terms coming from n'^ — r] 
mixing are (2m = md + m„) 

^(2) ^ a/3 md-mu ^ Q^^Qg^ ^ Q QQQg3 ^ 

4 nis — m 

and = ef^ + e% with [20, 21] 

(4) ^ A{2Ak -A„+Ar, + 2M|-) - 2M^iAK - A^) - 10247r^A"(3£^ + Lg) 
^ 247r2F2(M2 - M2) 

= £(2) (0.27 ± 0.09) = (3 ± 1) • 10"''' , (20) 
f4) _ -9Ze^(M| + Ak) + 167r2((2M|- + M^)Q{ - GM^Q-) _ ^ 4 .21) 

18\/37r(M2 - M2) ^ ' ' ^ ' 
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where A = Af |. - M^, A, = Mf log {Mf / n"^) , Q\ = 'SK^ - 6K^ + 2K^ + 2K^, Ql = + Kl^, and 
Zem — (-^^± ~ -^^o)/2e^F^ « 0.8. Numerical estimates are taken from Ref. [13] (which is based on 
the quark-mass analysis of Ref. [22] as well as from general dimensional analysis arguments for the Ki 
counterterms). 

Adopting the parametrization introduced in the first section, the QED corrections S^^j^.j and Jem are 
moved into the A^jm corrections (as well as the small s^^m effects), and will thus be dropped from Egs. fTT)) . 
It should be clear though that not all QED effects are moved into A em , since the physical charged particle 
masses occur inside the loop functions. In this respect, the separation into long-distance QED corrections 
and "purely" strong parts is somewhat ambiguous. 

From the general expressions, the isospin-breaking corrections for the slopes and for the form-factors at 
the origin can easily be obtained. The former will be discussed later on. For the latter, setting = in 
Eqs.lini), one finds 

ff^" (0) = 0.9775 - \/3(0.963e(2) + £(4)) ^ (22a) 
f^^^° (0) 0.9773 + \/3(0.963£(') + e^^') , (22b) 
f^"''^ (0) = 0.9773 - 0.0250£(2) , (22c) 
1+^''* (0) = 0.9775 + 0.0257£(^) , (22d) 
where is set to its physical value in the loop functions (as prescribed by unitarity). Setting instead 

= (2M^o + 2Af^± + M^o ~ 2Ml± )/3 , (23) 

as prescribed at this order, shifts all the form-factors down by a common 0.0004. 

Inserting e(^),e(^) as given in Eqs. (|19l20[) in Eqs. ([22)) . very precise estimates can be obtained. Unfor- 
tunately, the theoretical error on the isospin-symmetric Oijf") correction is significantly larger than the 
experimental errors on the Ki^ data. Specifically, the isospin-symmetric C(p^) result of Ref. [23] is, using 
the numerical estimates for the local terms found in Ref. [24] (see also Refs. [25,26]): 

(0) I o(p«) = ^i-op- + = -0.001 ± 0.010 , (24) 

OA 2 

Sloops = 0.0146 ± 0.0064, Sct = — ^(^[2 + C;^) = -0.016 ± 0.008 . (25) 

One should also keep in mind that the lattice estimate [27] for /^'^ (0) lo'pT)*^ = —0.016, and thus (^p""^*^ w 
—0.031. Further, it is only with this larger 0{p^) correction that Vus, as extracted from Ki^ data, satisfies 
the CKM unitarity. This seems to indicate that the error on Sct is under-estimated. As said, we will 
circumvent these difficulties by considering only two ratios of form-factors, to which wc now turn. 

3.2 The ratio r 

Once the long-distance QED corrections have been factorized into the A^m's (see Eqs. (l8l9p ). the form-factors 
satisfy 

fK+^° ( 2) fK+-K+ ( 2\ 

This relation is valid up to and including O [p^ , p"^ e'^^\ p"^ a) terms, since the isospin-exact 0{p^) correction 
drops out in the ratio. The dominant corrections then come from 0{p^s^'^\p'^a). Further, being momentum- 
independent, it tells us that at this order, the slopes of the CC form-factors are identical, as well as those 
for the FCNC form-factors. Since the same relation holds for scalar form-factors, we can write 

\CC — \K+Tr" _ \K"ir+ \FCNC — \K + it+ _ \K"7r" /07\ 



5 



From the two equalities in Eq. (|26p . one can form the double ratio r which is exactly 1 at this order: 



All the momentum-dependences cancel to first order in e'-^^ (and no further expansion in is needed). 

This is a striking prediction of ChPT. It can be understood from the quark diagrams for each transition, 
where the u or d spectator quark plays no role at leading orders, except for tt" — 7y mixing and of course long- 
distance QED corrections moved into Aem- In that respect, notice that all the QED corrections originating 
from the meson masses, kept in the form-factors, do also cancel out completely. 

At Ci(/e(2)), though t:^ — rj mixing still cancels in r, one could start to feel the effect of the spectator 
quarks. Naively, only non-local isospin-breaking contributions at two-loops could generate corrections (for 
instance "sunrise" -type graphs). Still, there is a good possibility that these effects also cancel out in r, 
at least to a large extent. A full two-loop computation would be required to check this conjecture, and is 
beyond our scope. For now, we will be satisfied by looking only at the behavior of ©(p^e'-^-') local terms 
for the vector form-factor (thus computing only tree-level wave-function and vertex corrections from the Ci 
Lagrangian of Ref. [28], in which we renormalize the Ci by F~'^). We observe that while there is a local, 
momentum-dependent correction to Ea. (|26p . 



1C Ap(2) 

ro+ = 1 + V3(£(2) + £(4) + £(6)) + (A(4Ci2 - 6C35) + (2Ci2 + Cgs + C90)) + O (a) , (29) 

128A^(2Afi + M,^) ^ Ci4 + Ci7 , C18 ^ , ^20 + C31 + C32 + 2C33 ^ 



all these terms cancel out in the double ratio r, apart from a small O (aZem) term 

r = 1 - ^^7f"" (2AC12 - q'C,o) + Oiis^'^r) (31) 

where C12 ^ a few 10~^ [26]. This indicates that indeed, there are strong cancellations at play in the 
ratio r between isospin-breaking corrections induced by the spectator quark. Therefore, in our numerical 
applications, we will take 

r = 1.0000 ± 0.0002 , (32) 
to account for possible {e^"^^)^ effects, 0{p'^a) corrections or r esidual ©(p^e^^)) 

non-local contributions. 

3.3 The ratio vk 

The ratio rK is defined as 



TK = +,o^+ ) [ = (1.00027 ± 0.00008) + (0.051 ± 0.01) e^^^ , (33) 
./+ (0) 



and parametrizes the isospin-breaking due to the initial kaon only, i.e. it is not sensitive to tt^ — r/ mixing. 
Compared to r, it cannot be expressed in simple analytic form, and is defined only for — 0. 

The errors are estimated as twice the shift induced by varying the 77 mass between its physical and 
theoretical value. Eg. ([^5]) . Still, as discussed in the previous section, the sensitivity to the spectator quarks 
is negligible at ©(p^e*^^^), and the above errors probably underestimate the full Oijf'e'^'^^) corrections. As 
before, a full analysis goes beyond our scope, and we consider only local contributions which are 

8 A 

(?-K)cT-0(pee(2)) = "^^^^'^^^'^ ~ 647raZe,„-p2Ci2 • (34) 
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BR(%) 
1?^^ 40.563(74) 



^e3 



27.047(71) 
5.0758(290) 
3.3656(280) 
0.07046(91) 



0.15454(29) 
0.10209(31) 
0.15889(30) 
0.10505(32) 
0.15454(29) 



Tabic 1: Result of the best-fit of ISTRA [29], KLOE [30,31], KTeV [32] and NA48 [33] data for the 
branching-ratios, as performed in [8], and phase-space integrals for each mode. 



We discard the O (p^a) correction, which is about 10% of the strong one if Cu ~ C34. Interestingly, 
the combination of counterterms Sct is exactly the one occurring in the isospin limit at 0{p^) for the 
form-factor at the origin, Eq. ((25)) . but is now suppressed by an additional e*-^-* factor. Using the numerical 
estimate 6ct — —0.016 ± 0.016, with an inflated error to account for the discrepancy between lattice (and 
CKM unitarity) and model estimates as well as for the neglect of 0(p^e^^^) loop contributions, and with e^^) 
from Eq. (fT9|) . we get a very precise estimate for vk'- 

TK = (1.00027 ± 0.00008) + (0.12 ± 0.07) e^^^ = 1.0015 ± 0.0007 , (35) 

which should thus include the bulk of ©(p^e^^^) effects. 



4 Numerical analysis 

We start in the next subsection with the estimation of the FCNC form-factor slopes and rare K decay 
phase-space integrals. Then, in the following subsection, we estimate the FCNC form-factors at the origin 
and discuss the K13 experimental situation. These values are then used in the third subsection to get the 
rare K decay matrix-elements, and finally, the last subsection deals with the long-distance QED corrections 
for TT^vP. 



4.1 Slopes and phase-space integrals 

The form-factor slopes are determined from Dalitz plot analyses (corrected for long-distance QED 
corrections). To an excellent approximation, the K'^ and Kl slopes are identical (see Eg . (l27t ) . A best-fit 
analysis of ISTRA [29], KLOE [30,31], KTeV [32] and NA48 [33] data was performed recently in Ref. [8], 
with the result 

= (24.82 ± 1.10) • 10-3 , /i -0.95 0.32 \ 

= (1.64 ±0.44) • 10-3 , p= 1 -0.43 . (36) 

Ao (13.38 ± 1.19) • 10-3 , \ 1 / 

From these slopes, the K^^ phase-space integrals are known to about 1/3%, see Table [TJ It should be noted 
that about a third of these errors is due to the poor quality of the fit (i.e., to scale factors) [8], originating in 
the strong experimental disagreement between Aq measurements. This will hopefully disappear in the near 
future. 

Let us relate these slopes to those of the FCNC form- factors. From Ea. ((27| . we know that they are equal 
(up to the electromagnetic effects, already included in A^m)- From Eqs. (fT7| . one can immediately get the 
isospin-breaking correction relating the CC and FCNC form-factor slopes 

\FCNC 

+ = 0.9986 ± 0.0002 . (37) 

K 

This O (p^£^^^) correction is very small, actually too small to really account for higher order effects. Indeed, 
the linear slopes only arise at O (p''), where they are dominated by the Lg counterterm, as for the pion 
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form-factor from which Lg is fixed. Since Lg is dominated by vector-meson exchanges, it is well-known that 
large SU (3) corrections at C(p^) are then needed to account for mp ^ rriK'- In our case, we of course do 
not try to fix Lg from the pion vector form-factor, but simply take the slopes as extracted from Kis data. 
We thus remain with the ©(p^e*-^-*) corrections, which we estimate from the measured K*~^ — K*^ mass 
difference as 



\FCNC . , 

A^*^ m^. 0(892) 



0.990 . (38) 



Therefore, to estimate the phase-space integrals relevant for rare K decays, we will rescale both the linear 
and quadratic slopes by 

^^FCNC 

+ P = 0.990 ± 0.005 . (39) 

This theoretical uncertainty will be seen to have a smaller impact on the phase-space integrals than the 
experimental uncertainties on the slopes themselves. 

For the scalar form-factor slope, the experimental information is less precise. Further, it is not so clear 
which resonance should play the dominant role as the presence of the k pole may blur the picture. We 
therefore assign a slightly larger error to account for ©(p^e*^^^) effects and set 

\FCNC 

" ^ ^0.99 ±0.01. (40) 

Nevertheless, this slope only matters for Kj^ 7r'^/i+/i^, and even there, the impact of this additional 
uncertainty on the final result will be very limited. 

Using the experimental slopes extracted from Eq. (|36p . we find 

J+ = 0.15269 ± 0.00028 ± 0.00007 , (41a) 

J,5' = lY'^ = 0.16043 ± 0.00030 ± 0.00008 , (41b) 

= 0.03766 ± 0.00010 ± 0.00003 , (41c) 

= 0.08624 ± 0.00068 ± 0.00009 , (41d) 

where the first error comes from the experimental error on the slopes, and the second from the unknown 
0(p^£^^)) terms. The correlation between 1^ and 1^ turns out to be of about 0.1% and can be safely 
neglected in evaluating B {Kl — ^ Tr^/i"*"/!") j^^py. Currently, the above errors are dominated by the experi- 
mental error on A^ (except I'^ for which Aq dominates). As a result, a reduction of about 50% for both the 
X'_f^ and Aq experimental errors, leads to a ~ 50% reduction for all the phase-space integral errors. 

4.2 Form-factors at the origin 

The Ki^ experimental branching ratios are given in Table [l] as obtained in the best-fit analysis of Ref . [8] 
(together with the tttt, tttttt modes and the Ks,l and lifetimes, tl = 51.173 (200) ns, T5 = 0.08958 (5) ns 
and — 12.3840 (193) ns). From these values, and after removing the long-distance QED corrections [8, 15], 
one finds 



IK. X f'^ ^ (0) lexp = 0.21645 (41) , (42a) 



K^3-- 


\Vus 


X 


J + 


(0)1 


= 0.21638 (55) 




\Vus 


X 




(0)1 


0.21678 (69) 




\Vus 


X 


jK"7r+ 


(0)1 


= 0.21554(142) 




\Vus 


X 


fK+TT° 


(0)1 


= 0.22248 (73) 




\Vus 


X 


fK+n" 
J+ 


(0)1 


= 0.22314(106) 



IK. X - (0) lexp = 0.22269 (60) . (42b) 
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The most precise strategy to get the FCNC form-factors at the origin from these data is to use the ratios 
r and rn as 



fK°7r+ (Q) 12 

IK. X ff'-- (0) I = rr^ '.^r iU^olr^i = 0.2107 ± 0.0010 , (43) 



,A"V" 1 _ IK. x/^ (0)|-,p 

IK. X (0) lexp 



IK. X (0) I = IK. X /f ■ (0) lexp = 0.2168 ± 0.0004 . (44) 

In this way, ©(p^e*-^)) corrections are under good theoretical controL Further, since the error on {r) 
accounts only for ^ 10% (3%) and ~ 28% (0%) of the final error, there is much room for improvements on 
the experimental side. 

Discussion: The strategy of using r instead of ro+ [4] to get |K. x f^°^° (0) | is more precise since, as 
discussed previously, higher order effects are under better control for r (see Eq. ip^ and ([31]) ). On the other 
hand, the disadvantage of using r is that we directly rely on K^^ data to estimate isospin-breaking effects. 

In this respect, it should be noted that neutral K^^ data are in much better shape than the charged 
Kg'^ data. Indeed, all the new generation Kaon experiments (designed for — K° mixing) have published 
the whole set of information relevant for neutral modes, i.e. branching ratios, Kl lifetime and phase-space 
measurements. In this way, K. x fi^ ^ (0) and its slopes are known with very high accuracy. At the 
moment this situation is not shared by the charged modes, for which the analysis relies on much older 
data, and for which the treatment of radiative corrections is unclear. New preliminary measurements have 
been announced by NA48 and ISTRA-h for the ratios B{K^^)/B {K+ — > vr+TT ), while KLOE has announced 
preliminary results for the absolute branching ratios [31] 

B{K+) = 4.965 (52) %, B{K+^) = 3.233 (39) % . (45) 

All these preliminary values are included in the best-fit analysis of Ref . [8] (Table [T]) . Unfortunately, the old 
PDG average for B [K^ tt+tt") is rather suspicious (see the discussion in Ref. [8])), while KLOE data are 
not precise enough to really compete yet. As a result, the fit to if^ data of Ref. [8] has a rather bad x^- 
This situation should soon improve, as both if+ tt+tt" and JsT^ modes are currently under study. 

The above comment can be made more quantitative by looking at the evolution of the experimental 
determination of the ro+ ratio as new K^^ data are announced: 

rK+^° (Q\ I ( 1.0328 ± 0.0039 Average K+ data, before Kaon 2007 [8], 

(ro+)cxp = ''"'^ +0^+ = < 1-0288 ± 0.0034 Average K+ data, after Kaon 2007 [8], (46) 



IK. X /f " (0) lexp 

|K. X f+ (0) lexp 1^ 2.0140 ± 0.0046 KLOE preliminary K+ data alone [31]. 



Within 2(7, these determinations agree with ro+ estimated at O (p*e^^') from the ~ r] mixing parameters, 
Eqs.dMSI]), 

(ro+)th = 1 + \/3(e(^^ + e^"^') = 1.0238 ± 0.0022 , (47) 

where errors are added quadratically. The 0{p^e'^^^) corrections (see Ea. (l29)) ) are accounted for in the above 
error, assuming they scale as e^^^/e^'*^ ~ e^'^^/e'^' 309(0. Therefore, waiting for the K'^^ situation to 
settle, we will also present final estimates for the rare K decay branching ratios based on the theoretical 
ro+ , Eq. (j47|) , together with K^^ data only, corresponding to 

|K. X ff-'" (0) Ith = 7^^^ IK. X ff^* (0) lexp = 0.2118 ± 0.0004exp ± 0.0005th • (48) 

Note that in principle, if Eq. ([171) is assumed to hold, one should perform the best-fit average of if^g and K^^ 
data to get |K. x ff (0) |av. = 0.21669 (36), and then use this value to estimate both FCNC form-factors. 
As this has only a very small impact numerically, and since we argued that the recent neutral K^^ data are 
in much better shape than the older data, we prefer not to use this average. 

^Comparing the local 0{p^) terms 5cTt Eq. ll25l l. to those for e(^) as well as for the spectator quark effects, Ea. ll29l l. one can 
note the larger numerical coefficients for the latter. Also, note that the spectator quark effects in ro+ start only at 0{p^e^^^). 
For these reasons, the naive counting could have to be slightly amended. 
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('^0+)th 


(^0+)oxp 






I 




r 


Future ? 




0.5173 ±0.0025 


0.5173 ±0.0025 


19 


43 


21 


17 




±0.0023 




2.231 ±0.013 


2.209 ±0.017 




76 


12 


10 


2 


±0.012 


"I 


0.7832 ± 0.0044 


0.7755 ±0.0058 




76 


12 


10 


2 


±0.0044 




0.1838 ±0.0011 


0.1821 ±0.0014 




71 


18 


9 


2 


±0.0010 




0.4210 ±0.0040 


0.4169 ±0.0045 




54 


37 


7 


2 


±0.0029 



Table 2: Final results for the rare K decay rate coefficients, in units of 10~^*'(|Kis|/0.225)^, for the ro+ 
estimates of Eg . (147)) and ([46]) (current average), respectively. For the latter, the approximate breakdown of 
the errors (in %) is also indicated. The last column shows the possible improvements achievable by a 50% 
reduction in the /_(?^^'^" (0), A'_,_ and Aq experimental errors. 



4.3 Coefficients for the rare K decay branching ratios 

For the n coefficients entering the rare K branching ratios, Egs. pU)) and one can in principle avoid 

propagating the error due to the Kl^ lifetimes by using, instead of Eqs. fl^ . the averages: 

{tl\Vus X ff-^ (0) He^p = 0.23978 (64) • 10-^ , (49a) 
{t+\Vus X f^^^" (0) He^p = 0.6141 (32) • 10-9 . (49b) 

Importantly, we do not use directly the iff 3 branching ratios, as is usually done following Ref. [4], because 
the radiative correction factors A/'^jy,^ have to be removed first (see Eq.®). 
From these averages, the Kl decay coefficients can be expressed as 



256,»8m*«„, ' *■ V " IK.. X /f (0) 



I cxp 



'w 



5 ■ 4/1 l^"'' {Tl\Vus X /+ (0) )oxp , (51) 

2567r5 sm dw V ^0+ / 



and similarly for 1%'^. The second form gives the smallest errors, i.e. optimizes the use of the experimental 



information. On the other hand, doing the same for the K'^ decay would increase the error because ro+ has 
a larger impact than : 



Kit 



5 . 4. lK.r T+ (r^lK. X /f - (0) |e.p)2 1+ (52) 
25o7r^ sm Bw 

G\M\a{Mzf ,8 , ^ fK+^« ,r.^\i^ ( tkV + , . 

= or-g q ■ 4/1 l*^"^! (^+1*^^^ ^ U (0) I )oxp , (53) 

25671^^ sm Ow \^a+ / 

Numerically, the first form leads to smaller errors. 

Using in addition the values quoted in Eg. pTj) for a (Mz) and sin^ 9^r, we get the results given in Table 
[21 for the experimental or theoretical ro+ estimates discussed in the previous Sectior|^. The breakdown of the 
error (in percent) is given for the (ro+)j,^p-based estimates, and shows that the experimental errors (on the 
if + lifetime on the K^^ form-factors at the origin f{0)Ke3 ^^'^ '^^ ^^'^ phase-space integrals I) completely 
dominate, so that there is much room for improvements. In particular, the last column shows the errors one 
would get if the errors on \Vus x (0) |cxp (from K^^, currently the most limiting), X'j^ and Ao were all 

reduced by a factor of two. 



^It should be noted that if the experimental uncertainty on sin'^ dyy was to be included in the and ft^ coefficients, their 
errors would increase by 14% and 11%, respectively. Since the present work concerns isospin-breaking effects, these uncertainties 
are left aside. Further, as said in the first section, they should be dealt with, together with higher-order electroweak effects, at 
the level of the short-distance Wilson coefficients [11]. 
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Figure 1: The QED correction to n^vvij), in %, as a function of the maximum energy of the 

undetected photon. 



FinaUy, of special interest for the search of new physics, the ratio of the two neutrino modes can be 
predicted with good accuracy: 



1 A4^T+(r+|KsX/f^-°(0)n 
r2 M^oXO(ri|K. x/f''-+(0)P),,p 



cxp 



= 0.2359 ± 0.0015 . (54) 



Further, the error is currently dominated by (r-|_|V„s x /jp ^ (0) | )cxp and could easily be reduced by a 
factor of two in the near future. 

4.4 QED radiative corrections 

The final piece needed is the long-distance QED correction to n^vv. For that, we have to combine 

the virtual photon correction given in Eq. (jl7dp with the real photon emissions. These have to be computed 
at O(p^a), leading to the radiative decay rate (see Appendix C for details and explicit expressions): 

v(T^+ +-\ GlM%a{Mzf\vA^ a f^'~''-^\ ,3 . , ^ ^ ,,,, 

r [K^ TT^VV^) = , . 4^ / dz\%JBR (2:,r^,Smax) ■ (55) 



2567r5 sin"* 9 



w 



The function JBR{z,rTy, E^s.^^) accounts for photon emission with energies up to £'max- Altogether, the 
IR-finite long-distance QED correction is 

a 1 /•(i-'-')' 

A_BM {Emax) = - -TT dzXl [ Jbr [z, r^, Sniax) - Jem {z , r^)) . (56) 

•^IiT Jo 

As shown in Fig[Tl Aem (-^max) is of less than one percent for reasonable -Emax, while the fully inclusive 
correction is Aem = —0.15%, and is thus comparable to those of Ki^ in magnitude. Still, being free of 
counterterms at leading order, it is precisely determined. Looking at Table[21 it is of the order of the current 
uncertainty on 

5 Conclusion 

In this paper, we have presented a theoretical strategy to extract rare K decay matrix-elements from Kg^ 
data to a few parts per mil. In summary, from the measured K^^ branching ratios and form-factor slopes: 

Step 1: The Kes phase-space integrals are directly obtained from the measured slopes, see Eg. pSj) 
and Appendix A. 
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Step 2: The Ki^ form-factors at the origin (modulo Vus^ kept as a parametric uncertainty) are 

L S -\- 

obtained after removing the QED corrections and the phase-space dependences from the K^^ 
branching ratios, and are separately averaged for the charged and neutral Kp^ decays, see Eqs. (|42p 
and (|49|. 

The outcome of these first two steps is provided by experimentalists, as in Ref. [8], taking into account 
the correlations among slopes, branching ratios and QED corrections. Then: 

Step 3: The slopes are rescaled to account for isospin-breaking, Eas. (|39l40p . and rare K decay 
phase-space integrals are computed (again, taking care of the correlations among slopes), see 
Eqs. pG]) and Appendix A. 

Step 4: Using the averaged values for the neutral and charged Ki^ form-factors at the origin 
together with the ratios r and Eqs. ([32)) and ((35)) . the rare K decay form-factors at the origin 
or, even better, directly the coefficients k are obtained through Eqs. (|51|52[) . 

Using this strategy, exceptional control over the hadronic physics is possible thanks to the two ratios 
of form- factors r and tk, for which higher-order corrections happen to be very suppressed. Overall, this 
leads to a significant improvement over the leading-order analysis of Ref. [4]. In addition, for the first time, 
uncertainties due to higher-order effects in the chiral expansion were carefully studied, allowing us to carry 
a thorough error analysis. Finally, QED corrections are accounted for, both for K^y, and rare K decays. 

The final results for the coefficients entering the rare K decay rates, Eqs. (|10ll3p . are collected in Tabled) 
As explained in the text, there is at present an open question related to the isospin-breaking effects observed 
in the ratio of the K^^ form-factor over the if^'g one (see Eq. ipS]) ). The current best-fit K^^ branching ratios 
relies on a number of quite old measurements and has a poor [8]. This problem will hopefully disappear 
in the near future, as the modes are currently under experimental study. Therefore, waiting for the K'^^ 
situation to settle, we consider the values in the first column of Table [2] as our best estimates. 

Compared to the matrix-elements used in Ref. [1], based on the leading order analysis of Ref. [4], the 
error for k+ (k^) is reduced by a factor of about 7 (4), respectively. Using these new values, together with 
the Wilson coefficient Eq.([3]) and the CKM parameters Eq.®, we find (adding errors quadratically) 

B {Kl -K^w) = (2.49 ± 0.39) x 10"" {Xt : 25%, : 3%, CKM : 72%) , , 

B {K+ Ti+vv (7)) = (7.83 ± 0.82) x 10"" {Xt : 20%, P„,c : 31%, k+ : 3%, CKM : 46%) ^ ' 

where we have taken A^m (^-max ~ 20MeV) « —0.30%. As the breakdown of the error shows, the long- 
distance uncertainties due to the matrix-elements are now negligible compared to the other (non-parametric) 
theoretical uncertainties. Similarly, for the direct CP-violating contribution to Kl n'^t^i^, the (non- 
parametric) error is now completely dominated by short-distance uncertainties in the Wilson coefficients uta 
and Hrvi Eq.(l5]). However, for these modes, it is the uncertainty over the indirect CP-violating contribution 
which is at present the most limiting, and for that, better measurements of B{Ks ~^ tt^^^^^) would be 
necessary. The impact of our new values on the B {Kl — ^ 7r°£+£^) predictions is therefore limited and we 
refer to [3] for the numerics. 

In conclusion, the hadronic uncertainties due to the matrix-elements of dimension-six FCNC operators 
are now well below the short-distance uncertainties in if — > iri/D and B {Kl — *■ ''^^^~^^~)]:,cpy -'^^ addition, 
being dominated by Ki^ experimental errors, they should become even more accurate in the near future 
since it is reasonable to expect slope measurements and K^^ data to further improve. Hadronic uncertainties 
will thus not hamper our ability to test the Standard Model with very high precision. 
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A Phase-space integral coefficients 

Expanding the phase-space integrals in powers of the slopes, and carrying out the integration, they can be 
expressed as numerical polynomials 



J = Co + ciX'+ + C2 (A? + A'l) + cgAf + cA'+K + '^sAo + , 
with the coefficients for each of the nine phase-space integrals collected in the table below: 



(58) 





Co 


Cl 


C2 


C3 


C4 


C5 


C6 




0.14084 


0.4867 


0.672 


2.276 


2.318 








0.09085 


0.2938 


0.474 


1.784 


1.751 


0.2055 


0.360 


r+o 


0.14480 


0.5006 


0.693 


2.354 


2.392 






— ^ 


0.09346 


0.3023 


0.489 


1.846 


1.809 


0.2120 


0.373 


J+ 


0.13958 


0.4721 


0.638 


2.068 


2.153 






7-0 n-V.A 


0.14603 


0.5156 


0.729 


2.586 


2.572 








0.03092 


0.2299 


0.453 


2.084 


1.892 








0.07665 


0.0691 


0.150 


0.804 


0.681 


0.5526 


1.198 



B Loop functions 

Explicit representations for the loop hmctions in the momentum range < < {M^i — -M^j)^ can be 
written in terms of a; = /M'^i,rT^ = M^j/Mj^;, r-y = M^/Mx' and 



2^1,2 



= l-rlTx, (i = l + rl-x, K = - 2{1 + x)rl + {I - xf 



(59) 



as 



m{x,r^) 



M? 



+ 2{rl - (1 - xf)x Xl P-X^ , (1 - riy 



■logr-TT - 77^ log 



2a;2 '^13 + X^ 



-4 (1 + r^) + ar ( 12877^1.^ _ log — fl + - 



and (with x = q^/Mj.^ and = M^± /M^±) 



Jem {x,r„) = (2-|-21ogrT, - logr^) ( 1 + ^log ^^-^p^"^^ + -^C{x,r„) , 



(60) 



(61) 



C{x,r^) = log^ 



2 f Z l + Att 

~2x 



, 2 /'^i - A^ \ 2 ^^2 - A^ \ 2 f ^2 + Xtt 

log I ) + log — - log 



2x 



2x 



2x 



+ 21og('^Vogf^U21og^ 



Z\ + J \ 2A,r 



Att ~ Z2 

2X^ 



Z2 - X^ 

-4Li2 



log 



Z2 + A^ 
2A, 



zi - A^ 

Zl + A,r 



(62) 



The IR-finite remainder obeys C(0,r7r) = 0. Also, one can check that Jbm (0, 1) = 0, as required by 
vector-current conservation. 



C Bremsstrahlung and QED correction 

For the bremsstrahlung process if + (P) 7r+ {K) v (pi) v (^2) 7 (fc), the O (p^e) amplitude is given by 

u{PuhA^-^'')''iP-) ^ (63) 



j^n ^ GF2ea{Mz) ^ 



V2 27r sin^ 0w 



2K°'P^' 



M^ + 2K-k M^-2P-k 



13 



since the vector form-factor is /+ = 1 at this order. Keeping > 0, this amphtude has to be squared and 
integrated over the four-body phase-space, decomposed as 

d<^4{P;K,k,p,,p^) = —d<^>3{P:K,T,k)d^2{T;p,,p^) . (64) 
Ztt 

The two-body leptonic integral can be done immediately, while the two variables for the three-body integral 
are chosen as the photon and pion energies. Introducing the reduced variables x = T^/M"^, z = 2Et^/Mk 
and u) = 2E^/Mk, the differential rate is 

dr{x,u),z) _ G]rMl-a{Mzf\yv\^ a ^ ZZ2 Xl z + lj {zi - lj) {2 - z - uj) - 2xz 
dxdwdz 2567r^ sin^ tt \ w a;^ 

_^,^^-£-a^L^^2.,-.-a; + r^l > (65) 
{f3-z-u}-r^) 

where terms leading to 0{r-y) contributions to the rate have been neglected. The integration bounds are 

< a; < (1 - r^f , 2r^ < a; < 1 - (r^ x/i) ^ a{x,y) - bix,y) < z < o [x. i,) + h {x, y) , (66) 

^ ^ (2-t.)(/3 + r2-u;) ^ J - 4r^^ [z^ - u + r^f - ir^x 

« (x, y) = , n — S , b {x, y) = 




2(l + r2-w) ' ^ '^^ 2(l + r2-w) 

There is no need to cut the x integral since the differential rate with respect to x is, in principle, an observable 
and is therefore IR-safe. 

To account for the experimental constraints, the photon integral will be limited to -Emax, i-t!- Wmax = 
Since there will always be some x for which 1 -|- — (r^ + ^/x) < uJmax, the phase-space 
boundaries become quite involved. To simplify them, we make use of the fact that the probability amplitude 
becomes imaginary outside the physical phase-space and write 

I dx duO (wmax — w) = Re / dx du) , (67) 

Jo J^Ty Jo J2ry 

as can be explicitly checked numerically. When working to leading order in Wmaxi the spurious imaginary 
part cancels out and after integration over z and w, we find the differential rate: 

TiK-^.-.,,) = 9l^^S^4^^ r'^'\xXlJ.,ix,r.,E^..) , (68) 



25677^ sin 6w ^ Jo 

where 

JBn{x,r^,E^^) = -(2+-\og^^J [log-- - -\og[j^J - -- — 

As can be checked, combining this with virtual photon corrections, the IR-divcrgcncc indeed cancels out. 

For numerical applications, we use the following numerical interpolation for the full correction (virtual 
-|- real photons) 

JbR {x, r^, i^max) - JeM {x, r„) = jo (Wmax) (l + X ji (Wmax) + O (x^)) , (70) 

jo (a) = 0.595 + 0.965 loga - 1.730a + 0.201a^ + O.lUa^ , 
ji (a) = -2.129 -0.0711 loga -I- 0.496a -I- 5.977a2- 5.585a^ . 

It is valid to better than 1% for all allowed values of Wmax- When integrated over x, it gives an approximation 
for the total QED contribution to the inclusive rate to within 5%, more than sufficient since the QED 
correction itself is very small compared to isospin-breaking corrections. 
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